Reduced density matrix functional theory at finite temperature. III. 
Application to the electron gas: Correlation effects 
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Based on our derivation of finite temperature reduced density matrix functional theory [1] and 
the discussion of the performance of its first-order functional [2] this work presents several different 
correlation-energy functionals and applies them to the homogeneous electron gas. The zero tem- 
perature limits of the correlation-energy and the momentum distributions are investigated and the 
magnetic phase diagrams in collinear spin configuration are discussed. 
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INTRODUCTION 

In Part I of the present work [1] we have presented 
the theoretical foundations of finite temperature reduced 
density matrix functional theory (FT-RDMFT). FT- 
RDMFT and its zero temperature counterpart employ 
the one-reduced density matrix (1RDM) rather than the 
density as central variable and are therefore also suitable 
for the description of systems subject to nonlocal exter- 
nal potentials. Furthermore, the availability of the eigen- 
values and eigenstates of the 1RDM, which are many- 
particle objects in contrast to the Kohn-Sham energies 
and states, opens up the possibility to describe phenom- 
ena that are difficult to access by DFT [3-9] . One promi- 
nent example are Mott insulators [9] which also establish 
the need for a finite temperature version of RDMFT. 

In view of the extraordinary success of the LDA in the 
context of DFT, it seems desirable to pursue a similar 
approach in the context of FT-RDMFT. If one consid- 
ers only local external potentials then there exists a big 
variety of methods to calculate the equilibrium proper- 
ties of the homogeneous electron gas (HEG) at either 
zero [10-13] or finite temperatures [14-27]. However, for 
an LDA in FT-RDMFT one needs the equilibrium free 
energy of the HEG subject to arbitrary nonlocal exter- 
nal potentials. These could in principle be calculated by 
Monte-Carlo calculations which, to our knowledge, have 
not been carried out so far. Accordingly, for setting up 
an LDA in FT-RDMFT, one needs to work with func- 
tionals that are approximate even for the HEG. Part II 
of this series [2] has therefore investigated the first-order 
functional of FT-RDMFT when applied to the HEG and 
it was shown that it was capable of reproducing a quali- 
tatively acceptable phase diagram including collinear as 
well as spiral spin configurations. The resulting critical 
densities, however, where far too large which emphasizes 
the need for approximate correlation-energy functionals 
in FT-RDMFT. 

The main goal of this third part of our work is the 
investigation of correlation functionals in FT-RDMFT 



when applied to the HEG. This will be done as follows. 
Firstly we will review the most relevant concepts of FT- 
RDMFT, discussing the possiblity of the formulation of 
a functional analoguous to the LDA. Secondly we will 
use Monte-Carlo results to investigate the relation be- 
tween correlation induced by electrons of the same spin 
compared to correlation resulting from the interaction of 
electrons with different spins. These results will rule out 
the most common class of approximations in unpolar- 
ized RDMFT as candidates for an accurate description 
of spin polarized systems and they will be used to guide 
the development of improved functionals. Thirdly, we 
will resort to the zero temperature situation and pro- 
pose a functional which for the first time not only re- 
produces the correlation-energy of the HEG over a wide 
range of densities but also describes the momentum dis- 
tribution qualitatively correctly. Finally, we will use the 
perturbative methodology as presented in [1] to derive 
a truly temperature dependent correlation functional in 
FT-RDMFT. 



THEORETICAL FOUNDATIONS 

We will now briefly review the most important con- 
cepts of FT-RDMFT [1]. A general quantum mechanical 
ensemble is described by its statistical density operator 
(SDO) D, which is a weighted sum of projection opera- 
tors on the underlying Hilbert space: 



(1) 



The corresponding 1-reduced density matrix (1RDM) is 
defined as the spatial contraction of all but one dimen- 
sions of D 



7(r,r') =tr{M+(r')*(r)}, 



(2) 



where the density p is determined via p(r) = 7(1", r). As 
the 1RDM, by construction, is hermitian we can write it 
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in spectral representation as 



7 (r,r') = ]>>0*(r> 4 (r), 



(3) 



where the eigenvalues {rii} are usually called occupation 
numbers (ONs) and the eigenstates {<fii} are called nat- 
ural orbitals (NOs) [28]. 

We have shown in [1] that the grand potential of an 
arbitrary system with possibly nonlocal external poten- 
tial in grand canonical equilibrium (eq) is uniquely de- 
termined by the corresponding cq-lRDM. The grand po- 
tential 17 of a system with 1RDM 7 can therefore be 
written as a functional with a universal part $[y], which 
is independent of the external potential of the system, 
and a part f2 ex t[7] which just depends on this external 
influence. 

fi[7]=3fr]+*W7]- (4) 
The external functional fl ex t [7] is given by 

n ext [ 1 ] = [ drdr'(v ex t(r',r)-»S{r-r'))i{r,r'), (5) 



where the chemical potential /1 governs the coupling to 
the particle bath. The universal functional ^[7] on the 
other hand is defined as 



£[7] = min tr{£>(T + W- 1/pinD)}) (6) 

= n k [ 7 ] - l/pSoh] + + fix [7] + n c [7] (7) 

where /? = l/(ksT) describes the coupling to the en- 
ergy bath with k,B and T being Boltzmann's constant 
and the temperature respectively. The different func- 
tional for the kinetic contribution fife [7], the Hartree 
and exchange contributions fi h [7] , &x [7] and the nonin- 
teracting entropy part So [7] are explicit functionals of 
the 1RDM and are given as 

fi fe [7] = J dr Hm (-?Pj 7 (r,r') (8) 

b] = \j dvdv'w{Y', r) 7 (r', r') 7 (r, r) (9) 

n B [ 7 ] = ~\f drdr'u;(r',r)7(r,r')7(r',r) (10) 

Sob] = ~Y1 ( n * ln ( n + i 1 - n ') ln ( x - ( n ) 

i 

Accordingly, one only has to approximate the remain- 
ing correlation functional £l c [7] . As the nonlocal nature 
of the 1RDM allows an exact description of the kinetic 
energy the correlation functional only contains contribu- 
tions from the interaction and the entropy which is ex- 
pected to simplify the construction of approximate func- 
tionals. 

The construction of the Kohn-Sham (KS) system in 
FT-RDMFT was shown to be particulary simple. The 



NO form the eigenstates of the KS Hamiltonian and the 
cigcnencrgics arc found by an inversion of the Fcrmi- 
Dirac relation: 
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(12) 
(13) 



With these findings we have all necessary tools at hand 
to describe the thermodynamic variables of a quantum 
system in grand canonical equilibrium. There are, how- 
ever, physical problems which require the formulation in 
terms of a canonical ensemble, e.g. the question of the 
equilibrium magnetic moment of a solid. In this situation 
the number of electrons is fixed and the equilibrium is re- 
alized by the state of minimal free energy F = E — 1/ j3S 
rather than minimal grand potential. 

As already mentioned in Part I of this work [1] the 
Hohenberg-Kohn theorem for canonical ensembles fol- 
lows exactly the same steps as for the grand canonical 
ensemble with the only difference that the external po- 
tential is uniquely determined by the 1RDM only up to 
an additional constant. Just as in the grand canonical 
case the free energy functional F[y] can be seperated in 
a universal part J c [7] and an external part F ext [7] where 
the universal part is defined as in Eq. (6) but the SDOs 
are now only defined on the TV-particle Hilbert space. 
The external part is defined as in Eq. (5), but without 
the chemical potential contribution. If we now assume 
our system to be in the thermodynamic limit, the func- 
tional for the free energy is simply related to the grand 
potential functional via 



F[ 7 ]=fi[ 7 ]+/iiV[7]. 



(14) 



This is due to the fact that in the thermodynamic limit 
the energy and entropy are identical for both canonical 
and grand canonical ensembles and therefore the uni- 
versal functionals $ c [y] and #[7] coincide. The exter- 
nal functional describes a one-particle contribution and 
therefore only depends on the 1RDM, regardless of the 
ensemble under investigation. 

After these considerations we can find the eq-free en- 
ergy F eq of a certain phase by a minimization of F[j]. 



F eq = minF[7] 



(15) 



The domain of minimization is given by the set of all 
ensemble-N-representable 1RDMS, i.e. lRDMs whose 
ONs fulfill 



0<n,-<l, J2 n * = N > 



(16) 



with a fixed particle number N. Further details of the 
numerical treatment of FT-RDMFT for the HEG can be 
found in [2]. 
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We have shown in Part II of this work [2], that the 
method of FT-RDMFT constitutes an efficient tool for 
the description of the equilibrium-properties of the HEG 
at finite temperature in FT-HF approximation. Eventu- 
ally, however, we would like to be able to describe cor- 
relation in spatially nonuniform systems, a task we will 
deal with in the following. 

CORRELATION IN FT-RDMFT 

Motivated by the big success of the LDA and LSDA 
in DFT, the formulation of a local-lRDM-approximation 
(LRDMA) in FT-RDMFT seems desirable. At this point 
it becomes important to point out the conceptual differ- 
ences between RDMFT and DFT. 

Local 1RDM approximation (LRDMA) 

In DFT, one only has to calculate the correlation- 
energy of the HEG as a function of the constant density 
p. In RDMFT, on the other hand, one has to calculate it 
for all possible density matrices 7(r — r') that are com- 
patible with translational invariance, i.e. for all possible 
momentum distributions n(k). The correlation-energy 
per volume of the HEG in RDMFT therefore becomes 
a junctional, e^ EG [y('k, R)], of n(k) rather than just a 
function of p. Under the assumption that such a func- 
tional was accessible, the LRDMA of RDMFT is defined 
in the following way. 

E LRDMA bM] ;= J d 3 fefBG[7(kjR)]) (17) 

where E LRDMA is the approximate correlation-energy of 
a nonuniform system and 7(k, R) describes the Wigner 
transform of the 1RDM, 

7(k,R) = yd 3 s7(R + s/2,R-s/2)e lk - s . (18) 

This procedure is similar to the definition of the LDA 
in superconducting DFT [29, 30]. While there are many 
LDA constructions conceivable that correctly reduce to 
the homogeneous limit, Eq. (17) is the only defini- 
tion that correctly reproduces the correlation energy of 
a weakly inhomogeneous electron gas [30]. 

A parametrization of e^ EG [n(k)] has not been carried 
out so far and remains an important task for the future. 
An extension of the LRDMA fomalism to spin-dependent 
systems and systems at finite temperature is conceptu- 
ally straightforward. However, the explicit calculation of 
the temperature-dependent e^ EG [n(k),T], although in 
principle possible by means of path integral Monte- Carlo 
techniques [31], is rather involved, due to the fermionic 
sign problem [32]. 



An alternative approach for further theoretical devel- 
opment in RDMFT and FT-RDMFT therefore consists in 
the derivation of approximate correlation functionals for 
the HEG and their implementation in a (FT-)LRDMA. 
We point out that a functional in FT-RDMFT should 
not only reproduce an accurate eq-free energy but also 
has to yield a good momentum distribution. 

The main focus of this work is therefore the presen- 
tation and investigation of several choices of correlation 
functionals in the theoretical framework of FT-RDMFT. 

Correlation-energy from DFT 

One choice which immediately suggests itself is the uti- 
lization of correlation functionals from DFT, i.e. func- 
tionals which just depend on the density p, rather than 
the full 1RDM: 

HcM = ficM. ( 19 ) 

We would like to point out that this approach exhibits 
an important intrinsic flaw. This becomes clear, if one 
reviews the minimization procedure for the case of lo- 
cal external potentials and seperates the variation over 
lRDMs in a combined variation over densities and a vari- 
ation over the class of lRDMs which yield a certain den- 
sity: 

F e „ = mill F["f] = min min ^[7] (20) 
7 p 7->p 

= min (Q ext [p\ + fillip] + Mc[p] 
p V 

+ min (fi fc [ 7 ] - k B TSoh] + ^H) ) ■ 

(21) 

As we assumed the correlation contribution to be inde- 
pendent of the particular form of the 1RDM as long as 
the diagonal, i.e. the density, stays the same we can take 
it out of the minimization over lRDMs. We find that 
the minimization over the lRDMs now only contains ki- 
netic, exchange, and entropy contributions. If the density 
p refers to a solution of the finite temperature Hartree- 
Fock equations to some external potential then we have 
shown in Part I of this work [1] that the constrained mini- 
mization of this combination of functionals gives the cor- 
responding Hartree-Fock momentum distribution. The 
resulting 1RDM will in this sense be "uncorrclated" . The 
severe problem with this result can be understood by con- 
sidering the zero-temperature limit. Any choice of corre- 
lation functional in FT-RDMFT which just depends on 
the density will lead to a step function as the momentum 
distribution for low temperatures [46]. This knowledge 
of the nature of the failure of DFT-functionals allows us 
to, at least partially, take account of these shortcomings 
and we will be able to study the effect of two representa- 
tive DFT-functionals later on in our investigation of the 
magnetic phasediagram of the HEG in FT-RDMFT. 



4 



We have seen in [2] that the first-order functional in 
FT-RDMFT was capable of reproducing a qualitatively 
resonable phase diagram. Apart from the incorrect mo- 
mentum distributions the main problem of these results 
was that the quantum phase transition between unpolar- 
ized and polarized configurations is instantaneous at zero 
temperature which is not correct for a quantum phase 
transition and that they occur at densities which are far 
too big (r s = 5.56a.u. compared to r s = 50a. u. ~ 75a. u. 
[11, 19]). We concluded that the noninteracting entropy 
functional is capable of describing the most important 
finite temperature features of the magnetic phase dia- 
gram, including a weakening of the ferromagnetic phase 
and the reproduction of a critical temperature T c , close 
to the Fermi temperature. As a first step in the develop- 
ment of correlation functionals in FT-RDMFT, it there- 
fore seems desirable to construct a zero temperature cor- 
relation functional which yields a critical density close 
to the expected one and then investigate the tempera- 
ture dependence induced by the noninteracting entropy 
functional. But before we are going to investigate such 
functionals we will first focus on an important property 
of the correlation-energy of a spin-polarized HEG. 

Spin-channel seperability 

The energy differences between different magnetic 
phases of the HEG become very close for low densities. 
An error in the spin-dependence of the correlation-energy 
might therefore result in a big error in the estimate of 
the critical density. Therefore, even in Monte-Carlo cal- 
culations, the estimates for the critical Wigner-Seitz ra- 
dius r c , where a quantum phase transition between an 
unpolarized and totally polarized configuration occurs, 
vary considerably (e.g. r s = 50a.w.:[ll] compared to 
r s = 75a.u.:[19]). An approximate functional in RDMFT 
would need to exhibit a very high accuracy to be able 
to reproduce the quantum phase transition in the HEG, 
which seems to be a very hard problem to solve. However, 
most physical systems of interest exhibit much higher 
densities and therefore the accurate description of the low 
density quantum phase transition of the HEG is rather 
irrelevant for the application of an LDA-type functional 
to real- world systems. This serves as an argument that 
for the implementation in an LDA also a medium accu- 
racy in the description of the correlation-energy would 
suffice. In the following, we will elucidate, why the use 
of a certain, rather general class of functionals will have 
severe problems to achieve an even mediocre accuracy in 
describing spin-polarized systems. 

When trying to describe spin-polarized systems in 
collinear configuration, the kinetic, the external poten- 
tial, and the exchange part can be seperated into a sum 
of contributions from only spin-up and only spin-down 
NOs and ONs. To describe the xc functional for a spin- 



polarized system, a common approach in RDMFT is to 
make the same ansatz of spin-channel seperability 

E XD [>y] = E xc [ ln , 7U ] = E s xc [ ln ] + E% c [>y u ], (22) 

where 7-^ and 74,4. are the diagonal elements of the 1RDM 
7 in spin-space and £ ,-S [7] denotes the functional which 
just depends on one spin-component of the 1RDM. 

We claim that such an approximation is intrinsically 
incapable of describing both spin-polarized and spin- 
unpolarizcd configurations together. The underlying rea- 
son for this problem is that an additive functional might 
describe the correlation contributions, which arise from 
interactions of the electrons of the same spin (i.e. equal- 
spin-channel correlation), well but will not be able to 
describe the contribution coming from the correlation 
of electrons of different spin (i.e. opposite-spin-channcl 
correlation). In the following, we will elaborate on 
this problem by considering a spin-polarized HEG. We 
will again rely on the PWCA[12] parametrization of the 
polarization-dependent correlation energies. 

For a collinear spin configuration, the fundamental 
quantities are the spin-up density and the spin-down 
density nj_ and the energy can be written as E(n,f,ni). 
If the assumption of spin-cannel seperability was valid, 
the following two relations would hold. 

E(n t , n±) = E(n t , 0) + E(0, r H ) =: E p (23) 

E(n t , ni ) = X - (E(n t , n t ) + E(n h nj) =: E u (24) 

In Eq. (23) the partially polarized system is given as a 
sum of two fully polarized systems, whereas in Eq. (24) 
one constructs the partially polarized one out of two un- 
polarized systems, hence the notations E p and E u . 

We can now investigate, whether or not Eqs. (23) and 
(24) are valid for the case of a partially polarized HEG by 
calculating the differences between the exact, i.e. Monte- 
Carlo, results and the expected results E p and E u . 

A p = E - E p (25) 
A u = E-E U (26) 

The results are shown in Figures 1 and 2. 

As we can see in Figure 1, A p vanishes by construc- 
tion for a totally polarized system (£ = 1). Decreas- 
ing the polarization then leads to a decrease of A p , i.e. 
E < E p . As in E p the opposite-spin-channel correlation 
is neglected, we deduce that it has to be negative. Con- 
sidering the unpolarized case in Figure 2, we see that 
again by construction A u vanishes, now for the unpo- 
larized configuration (£ = 0) and then increases with 
increasing polarization. This can be understood by real- 
izing that E u basically double counts the opposite-spin- 
channel correlation contribution, leading to E > E u . 

In RDMFT, the only correlation contribution comes 
from the interaction alone. We therefore would like to 
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FIG. 1: Energy differences A p from Eq. (25) for different 
polarizations 
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FIG. 3: Relative energy deviations 8 P from Eq. (27) for dif- 
ferent polarizations. 
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FIG. 2: Energy differences A u from Eq. (26) for different 
polarizations 



know, how the oppositc-spin-channcl correlation contri- 
bution relates to this correlation contribution W c . We 
can use the PWCA parametrization of the correlation 
contributions to the energy to calculate W c = E c — T c 
and investigate the following two fractions. 



5 U 



W c 



(27) 
(28) 



The results are shown in Figures 3 and 4. We see that 
by assuming spin-channel sepcrability, one will yield er- 
rors for the correlation-energy of up to 40%. However, a 
remarkable feature of Figures 3 and 4 is that the relative 
deviations over the whole range of considered densities 
do not vary strongly. Apparently, both equal- as well as 
opposite-spin-channel correlation are affected in the same 
way by a change in the density, leading only to a small 
change in their fraction. 

After these preconsiderations we are now prepared to 
investigate the performance of zero temperature RDMFT 
when applied to the HEG. 
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FIG. 4: Relative energy deviations 8 U from Eq. (28) for dif- 
ferent polarizations 



Zero Temperature RDMFT 

Most energy functionals in RDMFT try to incorpo- 
rate correlation via a modification of the exchange func- 
tional Eq. (10). The resulting functionals are accordingly 
called exchange-correlation functionals and can generally 
be written as 



E XC [ 1 )=E X [ 1 ]+E c [i] 
E xc [i\ = —x^2f(ni,nj) 



(29) 



dxdx' w (x, x' ) <p* (x' ) 4>i (x) 4>* (x) <pj (x' ) . 

(30) 

Choosing f x (m,nj) = niUj reproduces the exchange- 
only functional, neglecting correlation completely. 

The first approximation to f(rii,nj) was done by 
Miiller in 1984 [33, 34] leading to f Mueller {m,n 3 ) = 
yJriiTij . This Miiller functional was able to correctly de- 
scribe the dissociation limit of several small dimers of 
open-shell atoms but it overestimates the correlation- 
energy quite considerably. Inveigled by the simplicity 
of the form of the xc approximation and the primary 
success of the Miiller functional, Gritsenko, Pernal and 
Baerends developed what is now known as the BBC1, 
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BBC2 and BBC3 functionals [35]. The key difference in 
their approach is the different treatment of orbitals which 
were occupied or unoccupied in a Hartrcc-Fock solution 
(termed strongly and weakly occupied in the following). 
In addition, BBC2 and BBC3 effectively mix parts of 
the exchange and Miiller functionals. A very similar ap- 
proach was taken by Piris et al. [5]. Their approach 
differs slightly in the distinction of strongly and weakly 
occupied orbitals, removes parts of the self interaction 
(PNOFO) and tries to incorporate particle-hole symmetry 
(PNOF). These more elaborate functionals BBC1/2/3 
and PNOF/0 arc capable of reproducing good dissoci- 
ation energies as well as correlation energies. 

Whereas the previous functionals are mainly derived 
from physical arguments, Marques and Lathiotakis [4] 
pursued a different way by proposing a two-parameter 
Pade form for f(rii,nj). These parameters are then op- 
timized to minimize the deviation of correlation energies 
when applied to the molecules of the G2 and G2-1 sets. 
The resulting functionals are called ML and ML-SIC, ei- 
ther including or excluding self interaction. The ML and 
ML-SIC functionals achieve an unprecedented precision, 
reaching the accuracy of second-order M0ller-Plesset per- 
turbation theory for the calculation of the correlation en- 
ergies. An overview of several of these functionals as well 
as the corresponding energies for many molecules can be 
found in Ref. [36]. 

Another rather empirical functional, which will become 
important in this work, was derived by Sharma et al. [9] . 
They realized that both the exchange-only functionals as 
well as the Miiller functional can be seen as instances of 
a more general functional, namely the a or Power func- 
tional, described by f a (rii,nj) = (n.inj) a . As the Miiller 
functional underestimates the correlation-energy, one ex- 
pects to improve the results for values of a between 0.5 
and 1. This assumption proved to be valid [8], leading to 
an accurate description of the dissociation energy curve 
of H2 (a discussion of the performance of several DFT 
functionals for this problem can be found in [37]). 

These functionals were mainly constructed to describe 
small finite systems. However, because the correlation 
functional should be general, they should also work in the 
case of extended systems like the HEG. We will show the 
performance of some of the previously mentioned func- 
tionals when applied to the HEG later on in Figure 6. 

The BBC2/3 and PNOFO/ 1 functionals show a be- 
haviour similar to the one shown by BBC1 and ML-SIC 
performs even worse than ML. From the functionals men- 
tioned so far, only the a functional is capable of describ- 
ing the correlation-energy over the whole range of densi- 
ties considerably accurately. 

To make some statements about whether or not a given 
functional reproduces an accurate momentum distribu- 
tion, we will rely on the results by Gori-Giorgi and Zi- 
esche [38] . Three important properties of the momentum 
distribution derived from these results are listed in the 



following. 

(1) Nonzero occupation of all high momentum states 
[39] which can be understood by the help of the 
electronic cusp condition [1, 40]. 

(2) A discontinuity and rather symmetrical behaviour 
of the momentum distribution at the Fermi level 
(as also suggested by Landau-Liquid theory). 

(3) Depletion of low momentum states. 

All but one functionals mentioned so far succeed in re- 
covering property (1). This is due to the big value of the 
derivative ^g^l for m -> 0. All but the ML/ML-SIC 
functionals show a divergence for this derivative which 
will then lead to a partial occupation of all states. The 
ML/ML-SIC functionals exhibit only a very big, but fi- 
nite, derivative. Therefore, states of very high momen- 
tum are not to be expected to show partial occupation. 
When it comes to the description of property (2), only 
the BBC and the PNOF functionals recover a discontinu- 
ity at the Fermi surface. However, the overall behaviour 
around the Fermi level is not symmetric or qualitatively 
resembling the Monte-Carlo results. Furthermore, the 
size of the discontinuity, created by the BBC functionals, 
shows a wrong behaviour when changing r s . Although 
expected to decrease with decreasing density, it increases 
[41]. Property (3), a qualitatively correct depletion of 
low momentum states, is not fulfilled by the investigated 
functionals. Only the BBC1, BBC2 and PNOF function- 
als show a small decrease of occupation, but not nearly as 
much as prevalent in the exact momentum distributions. 

We believe that the depiction of low momentum states 
is an important physical effect which will become increas- 
ingly important in the low density limit and should there- 
fore be recovered by an RDMFT functional. Motivated 
by this credo, we are going to design an appropriate func- 
tional in the following. 

BOW functional 

We want to construct an xc functional which is capable 
of reproducing both occupation of high momentum states 
as well as depletion of low momentum ones. For the time 
being we do not focus on the discontinuity at the Fermi 
level. 

Considering the occupation of high momentum states, 
we will let us be guided by the success of the a func- 
tional and include a term of (jiirij) a in our functional. 
To achieve a depletion of low momentum states we then 
require our functional to have a vanishing derivative for 
m = rij = 1. In this way it is possible to reduce 
the occupation number of fully occupied orbitals with- 
out changing the energy. As the exchange contribu- 
tion is negative, this excess charge can be used to lower 
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FIG. 7: Momentum distributions of the 3D-HEG for r s G 
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FIG. 5: Occupation number dependence function / for the 
BOW and a functionals for different values of the parameter 
a. The Hartree-Fock function is reproduced by both func- 
tionals with a = 1. 



the xc functional. 



the energy. The variational principle will therefore lead 
to a groundstate where the orbitals are never fully oc- 
cupied. A possible choice of functional with vanishing 
derivative for full occupation could then be f(ni,rij) = 
((riinj) a — an.irij) / '(1 — a). We found, however, that 
this choice underestimates the correlation-energy consid- 
erably. The partially occupied states are given too much 
influence on the energy. We therefore introduce a simple 
counter-term to decrease this effect. It incorporates the 
inverse of the a functional, leading to our final choice for 
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FIG. 6: Correlation-energy of the three-dimensional electron 
gas using various correlation functionals from RDMFT. The 
black line denotes Monte-Carlo results in the PWCA [12] 
parametrization. 



dxdx'w(x, x')4>* {x')4>ii[x)4>* (x)<pj (x') (31) 



f BOW (ni,n,j]a) = (n^rij)" — arii-rij +a 

-a(l-n t n 3 ) 1/a (32) 

As a neccessary property, the BOW-functional leads 
to the reproduction of the exchange-only functional for 
uncorrclated momentum distributions 



f BOW (ni,nj; a)\n inj =o 
f BOW (rii, rij; a)| n4=nj=1 



(33) 
(34) 



Furthermore, we recover the exchange-only functional by 
choosing a — 1, i.e. f BOW (n i} rij;l) = f x (ni,n,j). We 
can now use the parameter a to tune the influence of 
correlation. We show f BOW (n u rij;a) and the function 
from the a-functional f a (rii, rij; a) = (riinj) a for several 
values of a in Figure 5. As we can see, a decrease in 
a leads to a bow-like shape which lead to our choice of 
name. 

The performance of the BOW functional when applied 
to the HEG in 3 and 2 dimensions in comparison to other 
prominent RDMFT functionals is shown in Figures 6 and 
9. 

Furthermore, we show several momentum distributions 
as resulting from the BOW functional in Figures 7 and 
8. 

In summary, the BOW functional shows remarkable 
success in describing the correlation-energy of the HEG 
over a wide range of densities (0.1 < r s < 100), including 
the range of metallic densities (1 < r s < 6). It manages 
to qualitatively correctly describe the depletion of low 
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FIG. 8: Momentum distributions of the 3D-HEG for r s e 
{5, 10}a.u. from the parametrization by Gori-Giorgi and Zi- 
esche (GZ) [38] and from the BOW functional for a = 0.61. 
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FIG. 9: Correlation-energy of the two-dimensional electron 
gas. The black line denotes Monte-Carlo results in the Attac- 
calite [42] parametrization. 



momentum states, an exact property known from Monte- 
Carlo results. 

Most of the functionals (Miiller, 
BBCl/2/3,PNOF/0,ML,ML-SIG) currently on the 
market are spin-channel separable which from our find- 
ings in the previous section will lead to a qualitatively 
bad performance when applied to the spin polarized 
HEG. The a and BOW functionals, on the other hand, 
offer a simple way out of the seperability dilemma 
because they incorporate a parameter. Making this 
parameter polarization-dependent makes the functional 
inseparable and offers an easy solution to the problem, 
given that the functionals also describe the partially 
polarized systems accurately. We have investigated the 
a and BOW functionals for partially polarized systems 
and found a good agreement of the correlation energies 
for the respective optimal parameters which we show in 
Table I. However, although both the a as well as the 
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0.61 


0.3 
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0.57 


0.63 


0.6 


0.58 


0.63 


0.7 


0.59 


0.64 


0.8 


0.61 


0.65 


0.9 


0.63 


0.67 


1.0 


0.66 


0.69 



2D 


(. 


a 


BOW 


0.0 


0.63 


0.66 


0.1 


0.63 


0.66 


0.2 


0.63 


0.66 


0.3 


0.63 


0.67 


0.4 


0.64 


0.67 


0.5 


0.64 


0.68 


0.6 


0.65 


0.69 


0.7 


0.66 


0.71 


0.8 


0.67 


0.73 


0.9 


0.70 


0.76 


1.0 


0.74 


0.80 



TABLE I: Best parameters to describe the correlation-energy 
of a HEG in three and two dimensions in the range 1 < r s < 
100 for the a and BOW functionals. 



BOW functional with polarization-dependent parameter 
describe the correlation-energy qualitatively correctly, 
they fail to predict a quantum phase transition between 
unpolarized and polarized phases. This is due to the 
fact that a smaller coefficient in the low density limit 
in both functionals leads to a bigger contribution to 
the correlation-energy and therefore, by reproducing 
the energy for intermediate densities well, they favour 
the unpolarized phase for low density, where the phase 
transition should occur. 

This is a backlash because the quantum phase transi- 
tion is a physical property of big interest. We will now 
propose a rather phcnomenological way to recover the 
quantum phase transition at zero temperature. From the 
good agreement of the BOW functional results with the 
Monte-Carlo results for fully polarized configurations, we 
deduce that one can, at least to some extent, describe the 
equal-spin-channel correlation effects by employing the 
exchange integral. We therefore propose to use a similar 
approach to include the opposite-spin-channel contribu- 
tions additionally. Our expression for this "trans-channel 
interaction energy" (TIE) modification reads with an ex- 
plicit mentioning of the spin index 

E BOW-T IE[l] = _^ f BOW {n ^ n . a . a P )K{iJ) _ 
ija 

C Y1 (( n ^ n ji) a (! -n^n j ^) a J K(i,j). (35) 





a p 


a u 


c 


BOW-TIE 


0.70 


2.0 


0.19 



TABLE II: Optimal parameters for the TIE-version of the 
BOW xc functional as defined in Eq. (35). 
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K(i,j) represents the exchange integral, corresponding 
to the NOs (pi and <pj and a p stand for the best pa- 
rameter for the description of the fully polarized HEG. 
The second term in Eq. 35 vanishes for a spin-polarized 
system and contributes increasingly with decreasing po- 
larization. The coefficients a u and c are fitted to re- 
produce a critical density closer to the Monte-Carlo re- 
sult while maintaining the good overall accuracy of the 
correlation-energy for different spin polarizations. The 
resulting parameters are shown in Table II and lead to 
an instantaneous phase transition at a critical density of 
r c si 28a. u.. 

We would like to emphasize again that we did not de- 
duce this opposite-spin-channel contribution from higher 
principles but rather postulated it to create a model func- 
tional which reproduces the critical density of the HEG 
more accurately. With a different choice for the inter-spin 
channel correlation energy, i.e. one which favours par- 
tially polarized configurations more strongly, one might 
be able to get rid of the instantaneous transition between 
unpolarized and polarized phases and reproduce a qual- 
itatively correct continuous quantum phase transition. 

This concludes our investigation of zero-tempcrature 
RDMFT and we are now able to investigate the effect of 
temperature in the framework of FT-RDMFT 



Finite Temperature RDMFT 

DFT-LSDA correlation 

Firstly, we are going to employ the correlation-energy 
of the electron gas as parametrized in the PWCA ap- 
proximation. Remembering the conceptual difference be- 
tween RDMFT and DFT, i.e. the knowledge of the exact 
kinetic energy functional, one might be tempted to re- 
move the kinetic contribution from the PWCA-functional 
before employing it as a correlation-energy functional in 
RDMFT. This procedure, however, would be incorrect as 
will be explained in the following. 

Based on our previous discussion of the intrinsic prob- 
lems of the utilization of DFT functionals in RDMFT 
we know that the minimizing momentum distribution 
will be "uncorrelated" , i.e. will reproduce the Hartree- 
Fock solution which becomes a step function n°(k) in the 
T — > limit. If we removed the kinetic contribution the 
correlation-energy functional would be given by the in- 
teraction contribution W PWCA alone and the free energy 
at density p would read 



F T ^° n fc [n°(k)] 



n x [n°(k)] + W t 



PWCA 



(P) (36) 




Because the momentum distribution is a step function, 
this expression reproduces exactly the Monte-Carlo re- 
sults without the kinetic contribution. This poses a prob- 
lem because removing the kinetic contribution from a 
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T(K) 

FIG. 10: Collinear phase diagram of the HEG for the DFT- 
LSDA functional. FM: ferromagnetic phase, PM: paramag- 
netic phase, PP: partially polarized phase. The dotted blue 
line denotes the Fermi temperature Tf- 



DFT-correlation functional yields an overall bad approx- 
imation. A straightforward remedy of this problem is the 
re-inclusion of the kinetic contribution to accomodate for 
the "uncorrelated" nature of the momentum distribution. 
Cashing in on our previous discussion we are therefore 
able to construct a FT-RDMFT functional which repro- 
duces the exact Monte-Carlo results in the T — > limit. 
The resulting phase diagram is shown in Figure 10. 

By construction, we recover the Monte-Carlo results in 
the PWCA parametrization for T — > 0. This leads to the 
reproduction of the continuous quantum phase transition 
at the correct critical Wigner-Seitz radius of r c w 70a.w.. 
In the radius range of 20a. u. to 70a. u. with an increase of 
temperature we encounter a continuous phase transition 
to a partially polarized state. A further increase of tem- 
perature then leads to a continuous phase transition back 
to the paramagnetic configuration. Considering the den- 
sity range defined by 70a. u. < r s < 90a. u., the situation 
becomes more complicated. Starting from the polarized 
configuration at zero temperature we encounter first a 
phase transition to a paramagnetic state, then another 
phase transition to a partially polarized state, and then 
a final transition back to the paramagnetic configuration. 
We attribute the appearance of the big partially polar- 
ized phase for temperatures above Tp to the fact that the 
1RDM is "uncorrelated" and therefore the noninter act- 
ing entropy functional underestimates the real entropy. 
It would be interesting to investigate, how a FT-LSDA 
from DFT would perform in the FT-RDMFT framework, 
but as mentioned before, the necessary Monte-Carlo re- 
sults do not exist so far. 

We will now turn to the investigation of a true RDMFT 
correlation functional. 
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FIG. 11: Collinear phase diagram of the HEG for the BOW- 
TIE functional. FM: ferromagnetic phase, PM: paramagnetic 
phase, PP: partially polarized phase. The thick dashed red 
line denotes an instantaneous phase transitions. The dotted 
blue line denotes the Fermi temperature Tf- 



BOW-TIE correlation 

Following from our considerations regarding zero tem- 
perature RDMFT we employed the BOW-TIE functional 
as defined in Eq (35) with the coefficients from Table II. 
The resulting phase diagram is shown in Figure 11. The 
critical density of the instantaneous transition increases 
with increasing temperature and the critical temperature 
is T w 47>. 

Apparently, the temperature independent BOW-TIE 
functional underestimates the effect of correlation as in- 
duced by temperature. To remedy this problem we will 
therefore in the following consider explicitly temperature- 
dependent functionals. We will again start with a DFT 
functional, namely the FT-RPA functional and will then 
propose a temperature dependent true FT-RDMFT cor- 
relation functional. 



FT-DFT-RPA correlation 



The correlation contribution to the grand potential in 
the FT-RPA is defined as 



{In (1 - W(q) X (q, !/„)) + VF(q) X (q, v a )} , (37) 
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T{K) 



FIG. 12: Collinear phase diagram of the HEG for the FT- 
DFT-RPA functional. The thick dashed red line denotes first- 
order phase transitions. Tf is too big to fit in the graph. 



where the frequencies v a are given by v a 
polarization propagator is defined as 



and the 



X(l,Va) = Xtfa'^a) +X4.(q,^a) (38) 

d 3 k 7i CT (k + q) - ?i CT (k) 



[ d A k ■ 
J v-ny lv a 



Mk + q)- £ff (k)) 



(39) 



In the case of DFT, the momentum distribution n(k) 
is given by the noninteracting Fcrmi-Dirac distribution. 
Results for the unpolarized HEG can be found in Refs. 
[22, 24]. We calculated the RPA correlation grand poten- 
tial for arbitrary polarization and temperature and show 
the resulting phase diagram in Figure 12. As the 1RDM 
is " uncorrelated" , at zero temperature the momentum 
distribution will be the same as in the FT-MBPT treat- 
ment, namely a step function. Therefore, we recover the 
critical Wigner-Seitz radius from zero-temperature RPA 
of r c sa 18a. u.. Compared to our general expectations 
the phase diagram only fails to correctly predict the crit- 
ical temperature T c is an order of magnitude smaller than 
the Fermi temperature Tp as 1300K of the corresponding 
critical radius as 21a. u.. We therefore deduce that the 
FT-DFT-RPA functional employed in the FT-RDMFT 
framweork manages to increase the critical Wigner-Seitz 
radius but overestimates the effect of temperature. 

The three functionals investigated so far were either 
intrinsically temperature independent or did just de- 
pend on the density alone. As a final approximation we 
will therefore now derive a correlation functional in FT- 
RDMFT which will be both explicitly temperature de- 
pendent and it will be a functional of the full momentum 
distribution. 
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FT-RDMFT-RPA functional 

We will utilize the perturbative expansion method as 
developed in Part I of this work [1]. We choose the subset 
of diagrams which refer to the RPA diagrams in FT- 
MBPT which are defined in Eqs. (37) - (39). Whereas 
the contributions Qk, Q ex t, So, and fl x only depend on the 
momentum distribution directly, fl r (and all higher-order 
diagrams which exhibit more independent momenta than 
noninteracting Green's functions) also depend explicitly 
on the Kohn-Sham energies. This becomes a problem in 
FT-RDMFT as we will show in the following. Using the 
fundamental one-to-one correspondence p f-> Vks, we see 
that p determines vks which determines e(k). Hence the 
e(k) are functionals of p and f2 r can be viewed as an 
implicit functional of e{k). To get a rough idea of the 
behaviour of this functional we parametrize the Kohn- 
Sham energies in the following way. 



e(m; k) 



2m' 



(40) 



where we choose m to be a variational parameter which 
can be interpreted as an effective mass. It will now be 
favourable to work with a reduced temperature t = T/Tp 
because it was shown [43] that for a noninteracting sys- 
tem the fugacity a = f3p only depends on t alone. From 
Eq. (13) we see that a dispersion relation as in Eq. 
(40) leads to a simple dependence of the fugacity on to, 
namely a(t;m) = a(mt), and the momentum distribu- 
tion becomes 



n(t, m; k) 



1 



1 



n(mt, 1; k). 



(41) 



The general effect of m on the momentum distribution is 
then that a decrease of m leads to a smoother momentum 
distribution. We are therefore able to investigate quali- 
tatively the behaviour of FT-RDMFT functionals under 
a change in n(k). Eq. (41) tells us that all functionals 
which just depend on the momentum distribution show 
the same simple dependence on m. 



VL k (t,m) 


= n k (mt,l) 


(42) 


VL ext {t,m) 


= Cl ext (mt, 1) 


(43) 


N(t,m) 


= N(mt,l) 


(44) 


S (t,m) 


= S (mt, 1) 


(45) 


n x (t, m) 


= fl x (mt, 1) 


(46) 


F(t,m) 


= Ek(mt, 1) + Q ext {mt, 1)— 






l/f3S (mt,l) + n x (mt,l) 


(47) 



It has to be noted that because of the temperature pref- 
actor of So, F(t,m) ^ F(mt, 1). The polarization prop- 
agator in RPA, on the other hand, shows a different be- 
haviour because of the explicit dependence on s(m; k). 



Let us now consider the limit of t — > while keeping the 
product mt, and therefore the momentum distribution, 
fixed, i.e. t — > 0,m cx 1/t. All contributions up to first 
order in the interaction stay invariant under this trans- 
formation. The free energy will change, because of the 
prefactor 1//3 in front of So- However, under the assump- 
tion that we started from a finite entropy this free energy 
change will be finite. The polarization propagator on the 
other hand will diverge because it exhibits a prefactor of 
to which leads to a divergence of fl r . 

Taking the previous considerations into account, we 
come to the conclusion that a straightforward inclusion 
of subsets of higher-order diagrams in the methodology 
derived in [1] will most likely lead to ill-behaved func- 
tionals, yielding wrong energies and momentum distri- 
butions. It has to be noted that this does not disprove 
the validity of the perturbation expansion of Opy]. It 
only shows that the utilization of a subset of diagrams 
in a variational scheme carries the danger of leading to 
a variational collapse. This fact is to be attributed to 
the total freedom of choice for e(k) by the inclusion of 
nonlocal potentials, a fact also recently pointed out in 
the context of GW [44]. As a general recipe for avoid- 
ing the problem of variational collapse one might model 
the perturbative expressions by approximations using the 
momentum distributions alone. We are going to investi- 
gate several first steps in this direction in the following. 

As we have seen before, the divergence in fl r stems 
from the fact that the eigenenergies £er(k) appear 
cxlicitely in the definition of \a ■ A first guess to remedy 
this problem is to fix these eigenenergies to the noninter- 
acting values. 



r (k)= e °(k) = ^- 



(49) 



We have implemented this correlation functional and 
found that it yields qualitatively incorrect results. The 
correlation functional in this approximation favours less 
washed out momentum distributions. 

This drawback lets us formulate a functional resting on 
exact properties of the polarization propagator. As a first 
step in our approximation process we will model x(Qj v a) 
as being frequency independent, much in the spirit of the 
COHSEX approximation [45]. Using the relation 



cm(k) 
de{k) 



-/3n(k)(l-n(k)), 



(50) 



the q — > limit of the kernel of x(q, 0) can be calculated 
up to second order to give 

n(k + q) - n(k) dn(k) 



e(k) de(k) 
0(n(k) - 5) V, .,'/-'/, t 



3e(k) 3e(k) 



Xa{t,m;q,v a ) = mx<r{mt, 1; q, v a ) 



(48) 



E;9 



dki 



(51) 
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We assume the ONs and eigenenergies to be point sym- 
metrical around the origin, i.e. n(k) = n(— k) and 
e(k) = e(— k). Accordingly, the second term in Eq. (51) 
vanishes in the integration over k. The polarization prop- 
agator in the limit q — > therefore becomes 



X.(q) ~P 



d 3 k 



n CT (k)(l- ?V (k)) + 0(g 2 ) (52) 



We now need to describe the momentum dependence 
for big q. Starting from the main expression for the RPA 
correlation grand potential Eq. (37) we use the spatial 
isotropy of the interacion and write 

fl r = ^(2tt)- 2 J dqq 2 {In (1 - W(q) X (q)) + W(q) X (q)} . 

(53) 

It is now possible to estimate the behaviour of £l r . If 
W(q)x(*l) is big, the logarithm in the kernel of ft r be- 
comes negligible. The main object defining whether or 
not W(q)x(q) is big or small is f3/(q 2 ). We therefore 
split the integral in Eq. (53) in two parts. One wherein 
q < q — \fj3 and one wherein q > q. 



r 

i 



(27T) 



r 
2 



dqq 2 (W(q) X (q)) 



4tt f vp 



roc 

-(2tt)- 2 / dqq 

P Jvp 



2 (W(q)x(q)) 2 



(54) 
(55) 

(56) 
(57) 



It can be deduced from Eq. (56) that in the limit T — > 0, 
i.e. (3 — > oo, Q° diverges unless %(q) behaves like 



X(q) — >oc q 
K < -2; 



(58) 
(59) 



We incorporate the findings from the investigation of 
the exact limits of \ m the following approximate k- 
functional 



X>(k)](q) = - 



P 



1 + g 2 /4 



^-p-n CT (k)(l -?v(k)) 



where 



(P 2 ) = 



d 3 k 
(2tt) 3 



(60) 
(61) 

(62) 



The prefactor of /3 ensures that Q r stays nonvanishing in 
the limit (3 — > oo. The occurrence of (p — (p 2 )) in this 
approximation shows that X £[ra(k)](q) will favour a more 
washed out momentum distribution. 
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FIG. 13: Correlation-energy of the three-dimensional electron 
gas for the ^-functional for various values of k. The black line 
denotes Monte-Carlo results in the PWCA [12] parametriza- 
tion. 



We calculated the correlation-energy in the zero- 
temperature limit for several values of k and show the 
results in Figure 13. The momentum distributions re- 
sulting from fi£[7] for k — 2.9 are shown in Figure 14. 

Unfortunately, the K-functional becomes spin-channel 
sepcrable in the T — ¥ limit. As in our discussion of 
zero temperature RDMFT we can make the parameter 




3.0 



FIG. 14: Momentum distributions of the HEG for r„ £ 
{1, 5, 10}a.tt. from the parametrization by Gori-Giorgi and 
Ziesche (GZ) [38] and from the K-functional for k = 2.9. 
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0.0 


0.1 


0.2 


0.3 


0.4 
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K 


2.9 


2.9 


2.9 


3.0 


3.0 


3.0 


t 


0.6 


0.7 


0.8 


0.9 
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3.1 


3.1 


3.1 


3.2 


3.2 





TABLE III: Best parameters to describe the correlation- 
energy of a HEG in three dimensions in the range 1 < r s < 
100 for the K-functional. 
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c 


k-TIE 


3.2 


0.24 



TABLE IV: Best parameters for the TIE-version of the k- 
functional as defined in Eq. (60). 
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FIG. 15: Collinear phase diagram of the HEG for the k-TIE 
functional. There are no first-order phase transitions. The 
dotted blue line denotes the Fermi temperature Tp. 



polarization dependent to describe the partially polarized 
HEG correctly (see Table III). However, for the same rea- 
sons as in the previous discussions the descreasing param- 
eter for decreasing polarization leads to a favourisation of 
the unpolarizcd phase in the low density limit and there- 
fore to the absence of a phase transition. We therefore 
propose a modification similar to the BOW-TIE func- 
tional as in Eq. (35) to recover a favourisation of the 
ferromagnetic phase for low densities. 



-TIE\ 



X ; "~[n(k)](q) 



d A k 



- x >(k)](q)- 

-n t (kK(k)(l-n t (kK(k)) (63) 



l + q 2 /k 2 F J (2tt) 

The parameters which lead to an accurate description 
of correlation energies for arbitrary polarization over the 
whole range of densities and a decrease in the critical 
density are shown in Table IV. 

We conclude our work by showing the the phase di- 
agram as resulting from the k-TIE functional with this 
choice of parameters in Figure 15. We see that the un- 
physical instantaneous quantum phase transition disap- 
pears and an increase in temperature slightly favours 
a partially polarized configuration over a paramagnetic 
one. In summary, the k-TIE functional constitutes an 
intrinsically temperature dependent true FT-RDMFT 
functional which describes the correlation-energy of a 
HEG at zero temperature with arbitrary polarization 
accurately, yields a qualtitativcly improved momentum 



distribution, as compared to previously used RDMFT 
functionals, and leads to a reasonable magnetic collinear 
phase diagram. 



CONCLUSIONS 

In this work we have investigated the possibility of 
the description of correlation effects in a HEG in grand 
canonical equilibrium in the newly introduced framework 
of FT-RDMFT. We have first focussed on the zero tem- 
perature situation and have shown that a certain class 
of functionals, namely spin-channel sepcrablc ones, will 
inevitably fail to describe the HEG for arbitrary polar- 
ization correctly. Using this knowledge we have then pro- 
posed a new functional, called BOW, which reproduces 
the correlation energies of a HEG at 3 or 2 dimensions for 
arbitrary polarization to unprecedented accuracy. It fur- 
thermore manages to describe the qualitative behaviour 
of the momentum distribution correctly. 

Wc have then turned to the description of the HEG at 
finite temperature and have shown that correlation func- 
tionals from DFT will lead to "uncorrelated" momen- 
tum distributions. Nevertheless we have utilized the zero 
temperature PWCA functional as well as the FT-RPA 
functional in FT-RDMFT and have investigated the cor- 
responding phase diagrams. 

Finally, we have investigated two truly momentum dis- 
tribution dependent functionals. Firstly the BOW func- 
tional and secondly a newly proposed K-functional which 
relies on exact limits of the RPA polarization propagator. 

This work being the first on the description of corre- 
lation in FT-RDMFT there are several possible proposi- 
tions for further research and improvements. Above all 
the need for Monte-Carlo calculations for the HEG sub- 
ject to nonlocal external potentials and/or finite temper- 
ature should be emphasized to get on the one hand a 
reliable measure for the accuracy of approximations and 
on the other hand the possibility of a real local-RDM- 
approximation. Further progress could be made by fo- 
cussing on the correct prediction of the discontinuity of 
the momentum distribution at the Fermi energy, possibly 
by the different treatment of states below and above the 
Fermi surface. 

In closing we would like to point out that the prob- 
lem of the description of different magnetic phases of the 
HEG lies at the core of the problem of an accurate de- 
scription of temperature dependent phases of real solids 
via a functional theory. However, the accurate descrip- 
tion of the low density regime in which the HEG under- 
goes a phase transition is complicated by the fact that the 
energy differences between distinct phases become very 
small over a wide range of densities. Real systems, on the 
other hand, show phase transitions at much higher den- 
sities where the energy differences change strongly with 
respect to the density. The BOW and K-functionals were 
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designed to reproduce the Monte-Carlo results for differ- 
ent phases over the full range of densities, including the 
physical ones. We therefore expect that they are suitable 
candidates for the description of correlation effects in real 
solids which shall be a task for the future. 
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